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Abstract In the absence of disorder, electrons can dis¬ 
play glassy behavior through supercooling the liquid 
state, avoiding the solidification into a charge ordered 
state. Such supercooled electron liquids are experimen¬ 
tally found in organic 0-MM' compounds. We present 
theoretical results that qualitatively capture the experi¬ 
mental findings. At intermediate temperatures, the con¬ 
ducting state crosses over into a weakly insulating pseu¬ 
dogap phase. The stripe order phase transition is first 
order, so that the liquid phase is metastable below Tg. 
In the supercooled liquid phase the resistivity increases 
further and the density of states at the Fermi level is 
suppressed, indicating kinetic arrest and the formation 
of a glassy state. Our results are obtained using classical 
Extended Dynamical Mean Field Theory. 


1 Introduction 

Glassy dynamics of electrons in the absence of disor¬ 
der [1,2] has long been elusive. Instead, quenched dis- 
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order seemed necessary for the kinetic arrest of elec- 
trons[3]. However, recent experiments on 9-MM' or¬ 
ganic compounds [4-6] found an exponential slowing down 
of the electron dynamics in a disorder-free environment. 
Inspired by these results, we proposed a model where 
the geometric frustration of the triangular lattice, aug¬ 
mented with long-range interactions, causes this glassy 
behavior. Using classical Monte Carlo simulations we 
showed that this model indeed reproduces the Arrhe¬ 
nius T ~ dynamics. [7] 

So how can glassy behavior arise without quenched 
disorder? The answer to this question can be found 
by considering ordinary structural glasses, like window 
glass. [1] A first order solidification transition is avoided 
by fast-cooling, leading to a supercooled liquid regime. 
Lowering the temperature of the supercooled liquid leads 
to the kinetic arrest and self-generated disorder. 

Similarly, in 0-type organics there is a solid phase of 
the electrons: the ground state is known to have stripe 
order. [4] The frustration caused by the triangular lat¬ 
tice and long-range interactions drive this transition 
first order, making it possible to form a supercooled 
electron liquid. Indeed, the cooling rate in 0-RbZn de¬ 
termines whether it enters the stripe ordered phase or 
the supercooled ’glassy’ phase. 

In these Proceedings we present theoretical results 
on the density of states and the resistivity of long-range 
interacting electrons on a triangular lattice, in the high- 
temperature liquid phase, the stripe solid phase and the 
supercooled liquid phase. We employ the methods of 
Extended Dynamical Mean Field Theory (EDMFT)[8- 
12]. Our results qualitatively capture the physics of the 
supercooled electron liquid, obtained by avoiding stripe 
order. 
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2 The model 


2.1 Long-range interactions on a triangular lattice 


The starting point is the model introduced in Ref. [7], of 
spinless interacting electrons on a triangular lattice. In 
the limit where the quantum hopping t is much smaller 
than the interaction strength V, the Hamiltonian is 
given by the interaction terms 


ij ^ / 



( 1 ) 


The electron density is fixed at one electron per two 
lattice sites. With only nearest neighbor interactions, 
this model corresponds to the triangular Ising model 
which is known to have a macroscopic ground state de¬ 
generacy [13]. We consider here instead the long-range 
Coulomb interaction between electrons, which lifts this 
degeneracy. 


Vij — 



( 2 ) 


where V is the nearest neighbor repulsion, and Vij is 
measured in units where the lattice constant is a = 1. 

The ground state of Eqn. (1) with Coulomb interac¬ 
tion is given by the stripe ordered state, with wavevec- 
tor k = M = ^0, ^^[7] At this M-point, the Fourier 
transform of the Coulomb interaction equals 


Ck=M = -1.40447C, 


( 3 ) 


from which one can hnd the stripe ground state energy, 
Eq = Vm/8. In the stripe ordered phase, the density of 
states consists of two delta peaks at ±|Vm|/2. Standard 
mean held considerations suggest a second-order phase 
transition at Tc = |Vm|/4. 


2.2 EMDFT 

A popular tool to study interacting electron systems 
is Dynamical Mean Field Theory (DMFT)[8], where 
the electron self-energy is found self-consistently as¬ 
suming it is momentum-independent. Extended DMFT 
(EDMFT) extends this idea to the electron polarization 
bubble. [9-12] In the classical limit t —>■ 0, where the 
polarization bubble becomes independent of frequency, 
the standard RPA expression for the polarization of an 
interacting system equals 

77k = 


(n - + I3{A + I 4 ) 
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Fig. 1 The self-consistent field /i as a function of tempera¬ 
ture T, as a solution of Eqn. (5). The density of states and 
conductivity follow from this result. In the low temperature 
limit, A cannot exceed the minimum value of the interaction 
in momentum space, VTnin- 


with A the only free parameter, n the average elec¬ 
tron density and 13 the inverse temperature. For classi¬ 
cal particles, the onsite density-density correlations are 
trivially TTioc = n — fi^. Because TTioc = 
arrive at the self-consistent classical EDMFT equation 

^ +/3(2i +Vk)' 


In the following, we will solve this equation for A as a, 
function of T, the results are shown in Fig. 1. 

Many other properties of the system, within the 
framework of EDMFT, can be derived from Z\. The 
local density of states p{uj) is given by the sum of two 
Gaussians,[9] 


p{uj) 


2 Z\ V 27 r V 



( 6 ) 


From the classical EDMFT results we can find the con¬ 
ductivity in a perturbative manner in leading order in 
the electron hopping t <^V. The dc conductivity gdc 
is then given by [9] 





pH^) 

AT cosh^ 


8 V Att 


( 7 ) 


The resulting density of states and dc conductivity are 
shown in Fig. 2 and Fig. 3, respectively. 


3 Pseudogap phase 

Given the solution of the EDMFT equations, we can 
analyse the different phases. At high temperatures, the 
system is a conducting electron liquid with resistivity 
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Fig. 2 The local density of states following Eqn. (6), in the 
correlated liquid phase. Around Tpc ^ 0.25V a pseudogap 
opens, suppressing the density of states at the Fermi level. 
This smoothly connects to the supercooled liquid phase where 
the density of states at the Fermi level is exponentially sup¬ 
pressed. 



Fig. 4 The free energy as a function of stripe order parameter 
m at a given inverse temperature 0 = IjT (at values 10, 12, 
14, 16, 18 and 20), following non-uniform EDMFT. From this 
we infer a first order transition temperature of Ts ~ O.OOF. 


4 The stripe ordered phase 


pdc 



Since the ground state is stripe ordered, one expects 
that at some finite temperature Tc the system enters a 
stripe ordered phase. To go beyond the standard mean 
field theory, we allow for the classical EDMFT equa¬ 
tions to break translational symmetry. To do so, we 
introduce in the electron density on a given site the 
stripe order parameter m such that 


(ni) = ^ + TO cos (Q • ri) = i ± m. 


( 8 ) 


The presence of stripe order reduces the on-site density 
correlations. 


Fig. 3 The dc resistivity following Eqn. (7). At Tmin ~ 0.14V 
the resistivity displays a minimum, indicating weakly insulat¬ 
ing behavior and a correlated liquid phase. Below Ts ~ 0.06F, 
the stripe order sets in with a first order transition (see Fig. 
4). In the stripe ordered phase, the resistivity is activated. 
This result is qualitatively similar to the experimental results 
on 6-RbZn, see Fig. lb of Ref. [6]. 


increasing with temperature, as can be seen in Fig. 

3. Upon lowering the temperature, however, around 
TpG « 0.25U the density of states opens a soft ’pseudogap’- 
like feature around the Fermi level, see Fig. 2. The 
pseudogap causes the system to become weakly insu¬ 
lating, with a minimum in the resistivity reached at 
T ■ « 0 14 

-^min ~ u.-L^T. 

The behavior of the resistivity is qualitatively simi¬ 
lar to experimental results in the organic compounds 9- 
MM'.[4-6] Both the resistivity minimum and the pseu¬ 
dogap signal the crossover to a strongly correlated elec¬ 
tron liquid. In this regime, seen in both experiments 
and our earlier Monte Carlo simulations, the system 
also develops short-range charge correlations. 


^\oC,i - (ll-Z 1^2 ) - ^ TO . 


(9) 


As a result, the EDMFT equation (5) is modified to 
1 2 1 


(I-to2)-i- b/3(Z\-hI4) 


( 10 ) 


where the self-consistent field A now depends on both 
TO and temperature T. The stripe order parameter to 
needs to be found from mean held theory, 

TO =-^ tanh Qto/3(Vq- b Z\(to,T))^ (11) 

which can be derived from the number density expecta¬ 
tion value. Thus for a given temperature the EDMFT 
equation gives Z\ as a function of to. Then, using the 
mean held condition Eqn. (11) we hnd the correspond¬ 
ing value of TO. 

Below T « O.IU we hnd two solutions for to: the 
symmetric liquid with to = 0 and an almost perfect 
stripe ordered phase with to « 0.5. Such a jump sug¬ 
gests a hrst order transition. To further corroborate this 
claim, we integrate the mean held equation to get the 
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free energy as a function of m, see Fig. 4. The free en¬ 
ergy profile clearly shows the presence of a hrst order 
transition, with the transition temperature Tg « 0.06t^. 

The transition into the stripe phase is so strong that 
immediately an almost perfect stripe order is created. 
With a density of states consisting of two delta func¬ 
tions, one expects the resistivity to be simply activated, 
Pdc ~ 4r cosh^ following Eqn. (7). The first order 
transition is therefore visible as a jump in the resistiv¬ 
ity, see Fig. 3, which is consistent with experimental 
results on organic compounds. [4-6] 

5 Supercooled liquid 

The presence of a first order transition, which is con¬ 
firmed within our classical EDMFT approach, allows 
for the possibility of a supercooled electron liquid with 
glassy dynamics. The symmetric m = 0 solution of the 
EDMFT equations remains metastable below Tc- Upon 
lowering the temperature, the pseudogap becomes stronger 
in the supercooled liquid which leads to an exponen¬ 
tial suppression of the density of states at the Fermi 
level, see Fig. 2. The resulting high resistivity (Fig. 3) 
is an indication of kinetic arrest: the motion of elec¬ 
tron basically stops and an amorphous glassy state is 
realized. Such an amorphous state is characterized by 
short-range correlations, as shown in our Monte Carlo 
results [7] and in X-ray diffuse scattering experiments [4]. 

When kinetic arrest itself arises, the system falls out 
of equilibrium and the manifestly equilibrium theories 
used here will no longer apply. At these lowest temper¬ 
atures we expect the self-induced disorder to cause a 
variable-range hopping type conductivity, yet a more 
detailed study of this regime remains a challenge for 
both experimental and theoretical future work. 

6 Conclusion and Outlook 

We have shown, using EDMFT methods, how kinetic 
arrest of electrons can arise in the presence of long range 
interactions on a lattice. Our work was motivated by ex¬ 
periments in 0-type organic compounds [4], though we 
expect our ideas to be valid in more materials. For ex¬ 
ample, many metal-insulator transitions are first order 
which allows supercooling. [14] In spin systems, recently 
the pyrochlore ’spin ice’ Dy 2 Ti 207 has been reinter¬ 
preted in terms of a supercooled liquid, [15] as well as 
the suggested spin-orbital liquid Ba 3 CuSb 209 [16]. In 
all these cases the combination of geometrical frustra¬ 
tion and long-range interactions seem to be the key to¬ 
wards understanding supercooled and glassy quantum 
liquids. 


Our results did not require any quenched disorder, 
which is the standard route for the generation of quan¬ 
tum glasses[3]. However, many features seem to be the 
same: the opening of a soft gap similar to the Coulomb 
gap[17,18], and the Arrhenius-type dynamics. [7] It is 
therefore an interesting question, to be addressed in fu¬ 
ture research, to what extent self-generated quantum 
glasses differ from quenched disorder-driven glasses. 
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